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Abstract 

The equation of state for tlie (3 stable nuclear matter is calculated numerically, and then the 
Tolman-Oppenheimer-Volkov(TOV) equation for the structure of the neutron star is solved in 
the fourth-order Runge-Kutta algorithm. It shows the mass and radius of the neutron star are 
functions of the central density of the neutron star and a maximum mass of 1.932 solar masses 
with a corresponding radius of 9.340km is obtained. Considering the equation of state of the 
nuclear matter must obey the causality, a new factor c is added in the nuclear potential energy 
formula. Therefore, with a new equation of state for the (5 stable nuclear matter when c = 0.15, 
a new maximum mass of 1.440 solar masses with a radius of 9.765km for the neutron star is 
obtained. Finally, the contribution of the cosmological constant to the structure of the neutron 
star is discussed, and we find the cosmological constant has minimal or negligible influence on the 
properties of the neutron star. 
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I. INTRODUCTION 



Since the first pulsar was discovered in 1967, it has been generally believed that the pulsar 
is a rapidly rotating neutron star endowed with a strong magnetic field. Thereafter, studies 
upon the properties of neutron stars gained wide concerns. 

Until now, many different views on the structure of the neutron star, especially the 
composition of the matter in the core of the neutron star, remain among theoretical physicists 
because it is difficult to ana^ 



n 



lyze the components of the neutron star matter from the present 
astronomical observations [l]. The core of the neutron star, where the density is very large 
and thence the kaon condensation could happen, possibly consists of the strange hadronic 
matter including hyperons and there could also be the large-size quark matter in it. In 
2006, it was concluded by Ozel that the condensates and unconfined quarks do not exist in 
the cores of neutron stars, because the mass and radius of the neutron star EXO 0748-676 
rule out all the soft equations of state of the neutron star matter [2]. However, Alford et al. 
stated that the quark matter can be as stiff as the nuclear matter, and the corresponding 
hybrid or quark stars can reach a mass of 2.0 solar masses^]. 

In the study of the structure of the neutron star, except the non-relativistic Brueckner 
theory Q, 3, the Quantum Had: 



great progress in the past 30 years | 
of state for the neutron star matter 
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3], which is also used to calculate the equation 



111 ]. To soften the equation of state for the neutron 



star matter, some research groups add self-interaction terms of mesons 121] in the Lagrangian 
density of QHD. At the meantime, considering the pair-correlation between nucleons, some 
scientists have established the relativistic Harteee-Bogliubov method 13|, [ij] . In addition. 



some others have developed the chiral hadronic model of QHD due to the chiral symmetry 



spontaneous breaking 
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16]. 



In this article, the nuclear model in Refs. [3] is utilized to calculate the equation of 
state for the neutron star matter and study further the properties of neutron stars although a 
similar work has been done by V. P. Psonis et al., where the relation between the asymmetric 
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enerery of the nuclear matter and the properties of the neutron star is discussed within 
the same model[19|. The equihbrium conditions that the chemical potentials of protons, 
neutrons and electrons in the neutron star matter will be discussed in Section |TT1 The 
equation of state for the neutron star matter and the requirement that the equation of state 
must meet to obey the causality will be discussed in Sections [1111 and HV] respectively, and 
the numerical results on the properties of neutron stars are analyzed in Section |Vl The 
probable influence of the cosmological constant on the neutron star will be studied in Section 
IVII The summary is given in Section IVIII It only needs some knowledge about quantum 
mechanics to calculate the equation of state with this model, so this work is suitable as a 
subject for the graduation thesis of an undergraduate physics major. Actually, this article 
is just based on Chang-Geng Liu's graduation thesis. 



II. (3 EQUILIBRIUM CONDITIONS 

Because a free neutron will undergo a beta weak decay process n ^ p + e~ + Ve, there 
are not only neutrons, but also a small fraction of protons and electrons in the neutron star 
matter. The neutron star is neutral. Hence, the density of the proton is equal to that of the 
electron, which means their Fermi momenta are the same as each other, kf p=kf^e. When 
the neutron decay process takes place, the electron capture reaction p + e~ n + Ve is going 
on. These two weak reactions will reach an equilibrium. This equilibrium can be expressed 
in terms of the chemical potentials for the three particle species 

A^n = A^p + fJ'e, (1) 

where the chemical potential for a particle species can be expressed as 

f^i = i^liC^ + m'^iC'^)K i = n,p,e (2) 

Thence, the Fermi momentum of the proton in the neutron star matter which reaches the 
beta equilibrium can be calculated as 

kf^J^ 3 
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III. THE EQUATION OF STATE FOR THE NEUTRON STAR MATTER 



Considering the influence of interactions between nucleons upon tlie equation of state, 
we mainly refer to Prakash's lectures on the asymmetric nuclear matter model to construct 



a simple nuclear potential energy formula 
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181]. The saturation density of the symmetric 



nuclear matter no = 0.16/m~^, and the binding energy per nucleon BE = E/A — 171^0^ = 
— IGMeV, where mjv represents the mass of the nucleon, and /cq is the nuclear compressibility 
that is applicable between 200 and 400MeV according to the experimental data and will be 
set to 360MeV in the discussion below. For the symmetric nuclear matter, the density of the 
neutron equals that of the proton, n„=np, with the total nucleon density n=n„ 
average energy per nucleon for the symmetric nuclear matter will be modeled as 17 



Hp. The 



n 5 zm^r 2 a + 1 

where u = n/riQ and a is dimensionless and A and B have the same dimension with the 

energy. The first term denotes the rest mass energy. The second term represents the average 
kinetic energy per nucleon with the corresponding Fermi momentum kf={3TT'^h^^)^ . At 
the saturation density, the average kinetic energy per nucleon is designated as < Ej > 
=22.1MeV , and then the second term can be rewritten as < Ej > u^. The last two terms 
in Eq. (jl]) indicate the mean nuclear potential energy per nucleon. From the properties of 
the saturation nuclear matter, three constraint equations can be obtained to determine A, 
B and a. 



!<B«>+^ + ^ = 0. (6) 

3 ^ 2 (T + 1 ' ^ ^ 

f <£»>+A + aB= ^, (7) 



which result in 



fco + 2 < E° > 
" 3<E'j> -WE' 

B = ^[i<E?>-fii?], (9) 
A = BE Ej> -B. (10) 
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If ko is determined, A,B and a are fixed accordingly. When ko = 360MeV, (T=1.922, 
A = -126.986Me1/ and B = -74.098Me1/. 
By the thermodynamic relation 



2d{e/n) de 



P = ''^^ = ''d^~'' ^^^^ 
the pressure in the symmetric nuclear matter can be expressed as 

p{n) = n,[\ < > + ^u" + (12) 

As for the asymmetric nuclear matter, the densities of the proton and neutron are not the 
same any longer. Let us represent the proton and neutron densities in terms of a symmetry 
factor a = """"^ with the proton density Up = ^-^n and the neutron density n„ = ^-^n. 

As the asymmetric part of the mean energy per nucleon is approximately proportional 
to the square of a, it can be written a.s\2^ 

E{n, a) = E{n, 0) + a^[{2^^ - 1) < > (ui - u) + Squ], (13) 

where Sq = 30MeV is the bulk symmetry energy parameter. Given the energy density 
e{n, a) = nouE{n, a), the pressure can be expressed as 

2 2 ^ 

p{n, a) = p{n, 0) + noa'^[{2'i - 1) < > {-ul - u^) + Squ^]. (14) 

Neglecting the electro-magnetic interaction between the electron and the nucleon, the 
energy density and pressure of the neutron star matter can be expressed as 

etot = £e + £{n,a), (15) 

and 

Ptot=Pe+p{n,a), (16) 

respectively, where Eg represents the contribution of the electron to the energy density and 
Pe to the pressure. 

eeikf)=eol (m^ + l)2M^(iM, (17) 

and 







Pe{kf ) = n—^ - Be = eo / {u^ + iy^u'^du (18) 

dn In 



With £0=^. 



a, n, kf^p and kj^e are all functions of fc/,„, so etot and ptot are both functions of /c/,„. 
Although the analytical expression between etot and ptot can not be deduced, the numerical 
form of the equation of state can be obtained. 

The fourth-order Runge-Kutta algorithm is employed to solve the TOV equation for the 
structure of the neutron star. 

dpjr) _ GM{r)e{r) p{r) ATrr^pjr) GM{r) 
dr c^r^ ^ £{ry M(r)c2 cV ^' ^ ^ 

and 

dM{r) 4nr'^e{r) 

where G is the Newtonian gravitation constant, e{r) is the energy density at distance r and 
p{r) the pressure, and M(r) is the mass of the matter included in the sphere with the radius r. 
To analyze the dependence of the mass and radius versus the central pressure, a sequence of 
central pressures are taken into the TOV equation. The energy density corresponding to any 
pressure is calculated with the linear interpolation method, and the computer programmes 
are written with Fortran. 



IV. AMENDMENT OF THE EQUATION OF STATE 

The square of the rate of the speed of sound in the nuclear matter over that of light 



satisfies 



18| 



(^)2 = = ^ = ^P/^^ ,21) 
c pc^ de de/du 

By analysis from Eqs. (jl]) and ( fT2i) . the speed of sound in the nuclear matter will exceed 
that of light at high nucleon density, i.e., the equation of state will violate the relativistic 
causality in this case. To make it obey the causality at the high nucleon density, the nuclear 



potential energy formula is revised as 



18| 



, . A B u , . 

rn,0 =-M + r, (22) 



where c is dimensionless and when it is set to zero, this formula returns to the unrevised 
case. Thereupon, the mean energy per nucleon for the symmetric matter turns into 

= mNC^+ <E^f>u^ + -U+ r, 23 

n 2 a + 1 1 + cu"-^ ^ ' 
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and the pressure changes as 



,2 „n B au''+^ 



,(n)=n,q<Ef>>ui + Y^-^^^L_^), (24) 

The parameters A, B, and a are fixed in the same way as in Section IIIIl but another 
parameter c has to be determined firstly. From our calculations, c is suitable between 0.1 
and 0.3. We set c=0.15 and obtain (t=2.246, A = -119AQ8MeV and B= - 80.562MeV for 
the nuclear compressibility ko = 360MeV. 

Fig. 1 shows the dependence of the energy density versus the nucleon density of the 
neutron star matter before and after amendment. At the low nucleon density, the curves 
fit together, while at the high nucleon density, they deviate from each other obviously. As 
can be seen in Fig. 1, the equation of state becomes softer after amendment. As is shown in 
Fig. 2, the curve of the equation of state after amendment only approaches the relativistic 
limit p=l£ and stays distinctly under the causal limit p=e, which indicates the equation of 
state after amendment obeys the causality at the high nucleon density. 
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Fig. 1 The energy density of the f3 stable nuclear matter versus the nucleon density 

before and after amendment. 
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Fig. 2 The equation of state for the (3 stable nuclear matter after amendment 

V. THE PROPERTIES OF NEUTRON STARS 

By solving the TOV equation, the masses and radii of neutron stars under different central 
densities are obtained. In Fig. 3, the dashed curve represents the dependence of the radius 
versus the central density of neutron stars when c = 0.15 and the solid one represents that 
when c = 0. Other than minimal deviation at the high nucleon density, the two curves fit 
together well basically, so the revised equation of state has little influence on the radius 
of the neutron star. However, the amendment of the equation of state does make great 
difference to the relation between the mass and the central density of neutron stars, as can 
be seen in Fig. 4, where the dashed curve is the case of c = 0.15, and the solid one is the 
case of c = 0. The maximum point of the solid curve corresponds to a maximum mass 
M^ax — 1-932M5„„ with a radius R — 9.370A;m, where Mgun is the solar mass, while the one 
of the dashed curve corresponds to Mmax—^-'^^^Msun and R— 9.765km. 
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Fig. 3 The radius versus the central density of neutron stars before and after amendment. 
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Fig. 4 The mass versus the central density of neutron stars before and after amendment. 



VI. THE INFLUENCE OF THE COSMOLOGICAL CONSTANT UPON 
THE PROPERTIES OF THE NEUTRON STAR 

If the cosmological constant A is positive, it has a upper limit of 4 x lO^^^s"^, and if 



negative, it has a lower limit of —2 x 10 ^^s ^ 2l|]. Considering the possible influence of A 

M{t)c^ 

M(r)c'^ 2GM{r)i- 

We take a series of values of A between the upper and lower limits into the revised TOV 
equation and solve it using the equation of state before the causality amendment. The 



on the structure of the neutron star, the factor [1 + ^^7:^73^] in the TOV equation in Eq. f[T^ 
should be changed into [1 + '^'"^ ^^^^ 



maximum masses in the unit of the solar mass and the corresponding radii in the unit of 
the kilometer under different values of A remain unchanged until the third place behind the 
decimal point. It implies the cosmological constant A has only minimal or even negligible 
influence on the properties of the neutron star. 

VII. SUMMARY 

The equation of state for the neutron star matter is calculated numerically, with the (5 
equilibrium conditions considered and the contribution of the electron included, and then the 
TOV equation for the structure of the neutron star is solved. It shows the mass and radius 
of the neutron star are functions of the central density. Considering the equation of state of 
the nuclear matter must obey the causality, a new factor c is added in the nuclear potential 
energy formula. Therefore, with a new equation of state of the neutron star matter when 
c = 0.15, a maximum mass and the corresponding radius of the neutron star are obtained. 
Finally, the contribution of the cosmological constant to the structure of the neutron star is 
discussed, and we find the cosmological constant has minimal or negligible influence on the 
properties of the neutron star. 

In the course of working on the graduation thesis on this topic, the student can learn the 
knowledge about the degenerate Fermi gas in quantum mechanics, know about the related 
contents of the special and general relativity and has an initial recognition of the structure 
of the neutron star. Besides, Chang-Geng Liu has learned to skillfully write computer 
programmes with Fortran to solve physical problems, utilize Origin software to draw figures 
and edit scientific articles with Latex. In addition, he has made great progress in writing 
and translating scientific articles. 
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